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| ■ Abstract 

We use the method of microscopic phase density to get general form of kinetic equation for the 

^— > | 

system of self-propelled particles with Vicsek-like alignment. The hydrodynamic equations are 

derived for the ordered phase and the ideal approximation is considered taking into account the 

mean-field force only. It is shown that the hydrodynamics of such ideal self-propelled fluid in main 

features demonstrates the dynamical transition from disordered initial state to the completely 

ordered motion. The problem of viscosity in hydrodynamic limit for the Vicsek-like model is 



discussed. 



INTRODUCTION 



The study of collective phenomena in systems of active agents is an important and fast 
developing field of statistical physics due to potential application to living matter dynamics 
[l|, |2|. This phenomenon is widely observed in nature at different levels of organization 
including the synchronicity in insects behavior, bird flocking, traffic and complex social 
behavior [l|, y-15] . These systems are essentially nonequilibrium and their dynamics is not 
strictly Hamiltonian due to information exchange which includes not only positions but also 
the velocities. 

One of the simplest models of the emergence of the coherent motion in active media is the 
(standard) Vicsek model (SVM) [6]. Basing on the simple rule it demonstrates the transition 
from disordered to the ordered state with varying the noise intensity. It has generated a series 
of publications |7H12| because of its simplicity and rich behavior. There are many questions 
that still need to be clarified for the model. One of the main open problem is the type of 
)hase transition. Still it is not known exactly either it is continuous [13l . ll4J or discontinuous 



15] . It is still the subject of ongoing debate even on the "experimental" level of numerical 



simulation. The mean-field approximations of network model 16] and of the Vicsek model 
itself 17| demonstrate the dependence of the type of the transition on the type of noise. 
Approaching this problem theoretically leads to the formulation of proper kinetics equation 
which adequately reflects the basic physical mechanism of self-organization and justifies the 
hydrodynamic equations proposed previously from phenomenological reasonings (3, Il8|, ll9[ 

Despite many results obtained in these approaches they lack the derivation of hydrody- 
namic equations which takes into account the specific features of SVM. The first attempt to 



do this was launched by Bertin et al. 
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20] and Ihle 2l|, |22] . In first works kinetic equa- 



tion was obtained within Boltzmann-like approach using two-particle collision integral. This 
gave the explicit expressions for the coefficients in hydrodynamic equations. But the results 
were quite different from numerical simulations, mainly because two-particle collisions are 
not sufficient to describe the model, and the usage of Boltzmann-like collision integral for 
active agents is not so evident, at least because the total momentum is not conserved. In 
the second works Liouville formalism was used. There the starting point was the equation 
for iV-particle distribution function with subsequent factorization to one-particle function. 
The obtained kinetic equation took into account multi-particle collisions but neglected cor- 



relations between particles, which are very important in the VM. The shear viscosity terms 
appeared in all approaches. But until now no experimental (numerical) evidence of the 
viscous shear was given in the VM. 

n 

In this paper we use well-known method of microscopic phase density [23] . In the frame- 
work of this approach the kinetic equation is constructed directly on the basis of equations of 
motion without reference to A-particle distribution function. Also we work with velocities 
v which form such additive quantities like flow density and pressure. We use the angular 
variables only to calculate final values, unlike previous studies where the angular variables 
were used directly and this led to highly nonlinear functions and increased the complexity 
of calculations. To the best of our knowledge firstly such approach for the self-propelling 



system with Vicsek-type of interaction was used by P. Degond in 24|, [25| . Although we also 



use the basic equation for the microscopic phase density functional as the starting point we 



proceed in physical approac 



in a standard way (see 



23 



i in order to get the collision terms and hydrodynamic equations 



[261) 

The paper is organized as follows. In Sec. [J we derive the equation for microscopic 
phase density functional which leads to the formal kinetic equation with the corresponding 
collision terms. Then we obtain the hydrodynamical equations. In Sec. [Til we consider 
the hydrodynamical limit of ideal self-propelled Vicsek-like fluid. The question about the 
existence of the shear viscosity for such fluid is the subject of Sec. IHIl There we study 
the viscosity of this fluid from the theoretical point of view using the equations derived in 
previous sections. We postulate the problem of numerical simulations of Coutte flow for 
Vicsek-like model and discuss preliminary results. The analysis of the results and problems 
for future studies are given in concluding section. 

I. THE EQUATION FOR THE MICROSCOPIC PHASE DENSITY FUNCTIONAL 

The standard tool for the derivation of the hydrodynamic as well as kinetic equations for 



dynamical systems is the method of microscopic phase density functional [231 ]: 

JV(z,t) = £<$(*-**(*)) (!) 

i 

where x = (r, v). It obeys the conservation law: 

d t N + vd r N + d v (irN) = 0. (2) 



The corresponding number density and the density flow are: 

pM(r,t) =fAf(x,t)dv, (3) 

j( m )(r,£)= fvA/-(M)dv, (4) 

where the superscript "m" stands for the "microscopic". We take the equation of motion 
for the Vicsek model [6j without noise term as following: 

dvi 



which descri 



[6 



dt -*•--" (5) 

Des the fact that the kinetic energy of a particle is conserved. The original Vicsek 



algorithm J6( is based on the alignment of the particle velocity to the average direction of 
the neighbors: 



W; 



Ej v j k (^ - r i, 



IEi v i^( r i- r i)l ' 

where K{vi — Yj) is the "microsopic" summation kernel and therefore: 



(6) 



u; i = 7v i xw i . (7) 

Usually, the Heaviside step function is used for this purpose with the characteristic region 
of several average distance between the particles. 

To avoid the difficulties in subsequent derivation connected with the non-additive struc- 
ture of Eq. (|6]) and to make the transition to the low density limit n — > trivial, we use 
another definition for w$: 

Wi = ^v i if(rj-r i ), (8) 

3 

i.e. the rate of alignment given by Eq. (jBJ) is proportional to the local flux of the nearest 
neighbors. Otherwise using of Eq. (|6]) leads to the essential multiparticle collision terms 



(see e.g. 2l|). These two options differ by the scalar factor 7. Therefore one can expect 
that despite the difference in microscopic equation of motion they lead in essential to the 
same macroscopic dynamics. Moreover, using Eq. (JS} makes possible the comparison of the 
results of this approach with the Boltzmann like kinetic approach of Bertin et al. [20]. 
Because of its additive structure Eq. (jSJ) can be expressed as the average of the MPDF: 

jM ( r , t) = ! dx' K(r - r') VM{x', t) = f K(r - r') j (m) (r', t) dr' . (9) 

Below we show that in such a case we get the equations similar to the usual case of molecular 
systems with pairwise interaction potential. 



The equation ([2]) takes the form: 
d t M{x,t)+vd T M{x,t)+^d v ( f dx' (vx v') x vK{r- r')N(x',t)N{x,t) j =0 . (10) 
From this equation it follows that the equation for one-particle distribution function 



Af(x,t) = nfi(x,t) is: 
d t fi(x,t) + v-d T f 1 (x,t) + -d v ( I K{t-v') (vx v') xwM{x,t)M{x>,t)dx'\ =0, (11) 



where Af(x, i)N{x', t) is connected with the pair distribution function f 2 : 



Af(x, t)N{x', t) = n 8(x - x') /i(x, t) + n 2 f 2 (x, x' , t) . (12) 

Using common definition of the correlation function g 2 = ji — /i /i we get: 



(8M8M) x ,x',t = Af(x,t)Af(x',t)-Af(x,t) Af(x',t) = n 2 g 2 ( y x,x',t)+n5(x-x') fi(x,t) . (13) 

Thus Eq. (TTTj) transforms into: 

(d t + vd r + d v F)f 1 (x,t) = I(x,t) , (14) 

where I(x,t) is the collision integral: 

I(x, t) — jndv I v x (v xV) K(r — r') g 2 (x,x',t)dx' ) . (15) 

Also we introduce the mean-field force F due to neighbors: 

F(x,t) =n7 / v x (v' x v) K{v - r') h{x',t) dx' = (vxj^)xv, (16) 

j K = j K(r-r')v' h{x',t)dx' = j ' K(r -r')j(r',t)dr' . (17) 

Using Eq. ( TT4l) the equations of motion for the basic physical quantities can be derived. 
The obvious conservation law for the number of particles is obtained from Eq. (I14p by 
integration over v: 

^ + divj = 0, (18) 

where 

p(r,t) =n f x (x,t)dv. 

The equation for the hydrodynamic velocity field u has the form: 

^^ + Ix~ 1 {pUlUl) = ~^ + :Fl + n S 5vJ{X)t)dW ' (19) 



where the pressure tensor P can be decomposed into scalar and traceless components in a 
standard way: 

*« =»/***» a* =*«,+«•«■ (20) 

Here 7Tjj is the traceless part which can be identified with the viscosity stress tensor 

iTij = n I 5vi5vj —Stj J fidv , (21) 

and 

Po = nJ i -^f 1 dv = ^(l-u 2 ) . (22) 

F(r,t)=njFf 1 (x,t)dv, (23) 

is the mean force density. In the hydrodynamic limit the spatial scale of the variations for 
/i is much bigger than that of the kernel K(r — r'). Therefore the microscopic kernel should 
be put to (^-function and j# = j. Thus the force can be decomposed as following: 

J 7 (r,t)=poJ-J.7r • (24) 



Also 



Note that in 24j the mean-field force was also obtained. It is proportional to the density 



gradient. But the representation in [24j did not provide clear distinction between this force 
and the dynamic term of the stress tensor because the authors used the field of unit vector 
^ — J/lJl °f the local flux. 

II. HYDRODYNAMIC LIMIT OF IDEAL SELF-PROPELLED VICSEK FLUID 

Neglecting all fluctuational terms in Eq. (1191) we get the equation of motion for an ideal 
self-propelling fluid as the continuous approximation of the Vicsek-like model: 

^_Zp + f(W)u (25) 

where po is given by Eq. (I2"2"j) . Although the first term in the right-hand side reminds the 
corresponding term of Euler equation for Newtonian ideal fluid. The second term as it follows 
from Eq. ( )24l) is the mean force due to neighbors. Note that this term has the structure of 
Landau theory like terms introduced phenomenologically in 7|. One can expect that taking 
into account the collision terms will renormalize the unit coefficients of the corresponding 




Figure 1. Density and velocity profiles with vq = — 1. 



power terms. The result Eq. ( )25|) re pre sents asymptotic form of the pressure terms in low 
noise limit. In the kinetic approach 20[ the corresponding term has the form (/i — £u 2 ) u 
where /z, £ are functions of density and the noise intensity. But the Boltzmann like approach 
of 20] is not applicable at low noise when u — > 1 since it is assumed that the density and 
the hydrodynamic velocity is low enough. 

Using the obtained equations for the ideal spp-ffuid it is natural to search for the running 
wave solutions of the form p(r — v£) , u(r — vt). Though Eq. fl5]) has sense only in dimen- 
sions more than 1 it is instructive to consider 1-dimensional case assuming the translational 
invariance along another directions. Here we give simple example of ID solution p(z) ,u(z) 
with z — x — Vq t. The equations (fT8l) and ( |25l) can be integrated explicitly: 

C 



P(z) 



u(z) - v 
where C is the constant of integration. The function u(z) is determined by the equation: 



(26) 



(1 - 2vl) In u(z) + v Q Oq - 1) ln(l - u(z)) + v (v + 1) ln(l + u{z)) = C z . 



(27) 



Nonsingular solutions with p(z) finite (see Eq. (|2"B"|) ) and < u(z) < 1 correspond to 
either Vq < or Vq > 1. The asymptotic behavior at t — > oo of Eq. ( 1271) depends on sign 
of Vq. If Vq < the order parameter field u — > (see Fig. [p. If v > 1 the transition to 
the completely ordered state u — > 1 takes place (see Fig. [2]). Note that the solutions with 




Figure 2. Density and velocity profiles with Vq = v2. 



u — > at t — > oo correspond to the situation where the disordered region has higher density 
than the disordered one. This can be interpreted as the disappearance of unstable state. 
This is expected result because of the consideration of the ideal regime of low noise where 
all collision terms are neglected. Note that these solutions are similar to the switching wave 
solutions in an active bistable medium [27 ]. 

In Fig. |3]we also show the solution of equations ( Tl8l) . ( 1251) under different initial states 
with periodic boundary conditions. It should be noted that these solutions asymptotically 
tend to the ordered state u(x,t) — > 1. 

Obtained results indicate that the ideal fluid approximation still retains the basic features 
of initial self-propelling system despite neglecting the collision terms. It shows how different 



terms introduced from the heuristic arguments earlier (see e.g. 



ay 



28|) arise. 



III. ON THE VISCOSITY OF SELF-PROPELLED FLUID 



In kinetic approach based on the Boltzman like arguments about the collision terms it 
is possible to get the equation of motion for the hydrodynamic fields [ill . |20| . Yet the 
appearance of the viscous terms seems quite unexpected at least if the Vicsek-like model is 
assumed as the basis for the derivation of the kinetic equations. At the microscopic level 
the standard Vicsek model |29[ is dissipativeless. Therefore one can expect the absence of 
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(a) u(x, 0) = 1 - cos 8 ( ff ) , p(x, 0) = 1 + cos 4 (ff ) 




(b) u(x, 0) = 1/4 , p(x, 0) = 1 + cos 4 (|f) 
Figure 3. Dynamics of the velocity u(x, t) and density p(x, t) under different initial conditions. 



viscosity at least in usual sense as the effect which decreases the momentum of the system. 

Accurate derivation of the collision terms along with the account for the noise terms will 

23 and 



be given in separate publication. Moreover different types of noise such as intrinsic 

the extrinsic one |9|] apparently have different influence on the order-disorder transition at 

least within mean-field approximation [14|, |l6|, [30 1 . 



Here we show that in the local equilibrium approximation the off-diagonal part of viscous 
tensor Eq. ( I21~j) vanishes. Indeed, as it follows from 17|, |24], |3JJ in the local equilibrium 
approximation for the ordered phase the distribution function can be approximated as fol- 
lowing: 



./ 



(o), 



Ce a 



(28) 



where C and a are some functions dependent on the local value u(r) of the order parameter. 
Obviously, it can serve as the interpolation between homogeneous distribution function 
/ = ^- in the disordered state and the formal series representation of the distribution 



function near transition point to the ordered regime [20]: 

f(v,e,t) = ^ (±- + -u(r,t) cos 9 + 

n \ 2ir 7i 



In such approximation the viscous tensor is: 

ViSvi-Sn^- f {0) dv (29) 

The off-diagonal components % ^ j vanish 



7T (0) 



tt 12 = tt 21 = j S Vl Sv 2 Ce acos0 d9 = C I "(cos 9 - u) sin 9e acos8 d9 = (30) 

according to the symmetry. The diagonal part of the tensor is determined by the quantity: 

n(°> = )- ![5vl - 8v 2 2 ]f {0) dv = ]- f[(cos9 - uf - sin 2 9]f^d9 . (31) 

which is not zero in general. Thus in this approximation the viscosity tensor is represented 
as: 

*$> = n<°> r _° 1 (32) 

Obtained results show that in this respect spp-ffuid behaves differently from molecular 
liquids. For the ordinary fluid the local equilibrium is given by the Maxwell distribution 
and ti\, = due to basic properties of the collision integral [23| . We have shown that the 
shear viscosity is absent in the Vicsek-like fluid in local equilibrium approximation given by 
Eq. f j28l) . Previous studies showed the viscosity terms in the hydrodynamic equations. So, 
this question needs to be resolved. The simplest way to check the existence of viscosity is 
a numerical simulation. To do this one can perform simulations of self-propelled particles 
system subjected to shear (as example [32]) and study the profile of the velocity field u that 
appears in such system. In the usual viscose liquid (Couette flow) linear profile (gradient) 
appears. Expected zero viscosity results in the absence of the nontrivial profile of the order 
parameter, except, probably, very thin layer close to the aligning boundary, caused by the 
specific mechanism of interaction. 



We have performed a first simple check for a small system 33J and preliminary results 
show the absence of the order parameter profile. But the proper numerical answer to this 
problem is a question of our future study. 
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IV. DISCUSSION 



In this paper we use the standard method of microscopic phase density to derive the 
general form of the kinetic equation for the Vicsek-like model. We obtain the form of the 
collision terms which in principle gives the possibility to obtain the closed kinetic equation 
for one-particle distribution function f\. Such an approach is more closely connected with 
the equation of motion and could be easily modified to include different kind of velocity 
constraints and noise perturbation to dynamics. We have considered the simplest mean-field 
approximation and checked that the corresponding hydrodynamic equations have solutions 
displaying the transition to the ordered state. Note that we considered the rule of alignment 
to the local flux of the neighbors, but not to the average velocity as in original Vicsek model. 
By this the calculations are simplified because of the additivity and collision integrals depend 
on pair correlation function only. Although this simplification seems to be appropriate from 
the physical point of view, the question of mathematical equivalence between these two ways 
remains open and proper numerical study is required. Usage of the original Vicsek rule of 
alignment to the average velocity in principle lead to essentially many-particle collision terms. 
Thus in such an approach it is possible to consider in unified terms the kinetic equation 
obtained in Boltzmann-like derivation scheme [ll|, 



|20( and master equation approach 21] 
which lead to essentially many-particle collision terms. 



Also we have considered important question about the viscosity of the Vicsek-like model. 
The calculation of the viscous tensor in the local approximation shows vanishing of the 
diagonal components of the tensor. We have checked this by looking at the velocity profile 
in a small test system. However proper numerical proofs of this result are required. The 
'orm of distribution function was taken as Boltzmann-like exponent from previous works 



17 



31] though more rigorously it should be evaluated self-consistently on the basis of the 



kinetic equation. These and other questions will be considered in separate work. 
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